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> ■ ABSTRACT 

o : 

, Aims. We present a qualitative analysis of key (but yet unappreciated) linear phenomena in stratified hydrodynamic 

' ■ Keplerian fiows: (i) the occurrence of a vortex mode, as a consequence of strato-rotational balance, with its transient 

7—i ' dynamics; (ii) the generation of spiral-density waves (also called inertia-gravity or gQ waves) by the vortex mode 

^SJ ' through linear mode coupling in shear flows. 

Methods. Non-modal analysis of linearized Boussinesq equations were written in the shearing sheet approximation of 
accretion disk fiows. 

Results. It is shown that the combined action of rotation and stratification introduces a new degree of freedom, vortex 
mode perturbation, which is in turn linearly coupled with the spiral-density waves. These two modes are jointly able 
^ , to extract energy from the background flow, and they govern the disk dynamics in the small-scale range. The transient 

5^ . behavior of these modes is determined by the non-normality of the Keplerian shear flow. Tightly leading vortex mode 

' perturbations undergo substantial transient growth, then, becoming trailing, inevitably generate trailing spiral-density 

' waves by linear mode coupling. This course of events - transient growth plus coupling - is particularly pronounced for 

I— I, perturbation harmonics with comparable azimuthal and vertical scales, and it renders the energy dynamics similar to 

■ the 3D unbounded plane Couette flow case. 

04 ' Conclusions. Our investigation strongly suggests that the so-called bypass concept of turbulence, which has been recently 

^ ' developed by the hydrodynamic community for spectrally stable shear flows, can also be applied to Keplerian disks. This 

OO , conjecture may be conflrmed by appropriate numerical simulations that take the vertical stratiflcation and consequent 

. mode coupling into account in the high Reynolds number regime. 

^ ■ 

^.f^ . Key words, accretion, accretion disks - turbulence - hydrodynamics 

d ; 

J — 1. Introduction totic flow stability to the study of short time behavior. It 

, was established that asymptotically /Rayleigh stable flows 

. According to classical fluid dynamics, unmagnetized disk allow for linear transient growth of vortex and/or wave 

^ I flows in Keplerian rotation (more generally: with angular mode perturbations (cf. Gustavsson 1991; Butler and Farrell 

^ I momentum increasing outward and with no extremum of 1992; Reddy and Henningson 1993; Trefethen et al. 1993). 

^ . vorticity) are spectrally stable; however, there is irrefutable This fact incited a number of fluid dynamicists to examine 

C3 ■ observational evidence that such disks have to be turbu- the possibility of a subcritical transition to turbulence, with 

lent. Due to this apparent contradiction, disk turbulence the linear stable flow finding a way to bypass the usual route 

is often considered as some sort of mystery. However, an to turbulence (via linear classical/exponential instability), 

analogous dilemma that existed in laboratory/engineering On closer examination, the perturbations reveal rich and 

flows has been solved by the hydrodynamic community complex behavior in the early transient phase, which leads 

in the 90s of the last century, where a breakthrough was to the expectation that they may become self-sustaining 

accomplished in the comprehension of turbulence in spec- when there is nonlinear positive feedback, 
trally/asymptotically stable shear flows (e.g. in the plane gascd on the interplay of linear transient growth and 

Couette flow). nonlinear positive feedback, a new concept emerged in 

Let us briefly recall the essence of that breakthrough the hydrodynamic community for the onset of turbulence 

(see Chagelishvili et al. 2003 for details) . Traditional stabil- in spectrally stable shear flows and was named bypass 

ity theory followed the approach of Rayleigh (1880) where transition (cf. Boberg & Brosa 1988; Butler & Farrell 

the instability is determined by the presence of exponen- 1992; Farrell & loannou 1993; Reddy & Henningson 1993 

tially growing modes that are solutions of the linearized Gebhardt & Grossmann 1994; Henningson & Reddy 1994 

dynamic equations. Only recently has one become aware Baggett et al. 1995; Grossmann 2000; Reshotko 2001 

that operators involved in the the modal analysis of plane Chagelishvili et al. 2002; Chapman 2002). The bypass sce- 

shear flows are not normal, hence that the corresponding nario differs fundamentally from the classical scenario of 

eigenfunctions are non-orthogonal and would strongly in- turbulence. In the classical model, exponentially growing 

terfere (Reddy et al. 1993). For this reason, the emphasis perturbations permanently supply energy to the turbulence 

was shifted in the 90s from the analysis of long time asymp- and they do not need any nonlinear feedback for their self- 
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sustenance, so the role of nonlinear interaction is just to 
reduce the scale of perturbations to that of viscous dissi- 
pation. In the bypass model, nonlinearity plays a key role. 
The nonlinear processes arc conservative, but in the case 
of positive feedback, they ensure the repopulation of per- 
turbations that are able to extract energy transiently from 
the mean flow. The self-sustenance of turbulence is then 
the result of a subtle and balanced interplay of linear tran- 
sient growth and nonlinear positive feedback. Consequently, 
thorough examination of the nonlinear interaction between 
perturbations is a problem of primary importance, and the 
flrst step is to search and to describe the linear perturbation 
modes that will participate in the nonlinear interactions. 

Such linear transient growth is also at work in rotating 
hydrodynamic disk flows; however, the Coriolis force causes 
a quantitative reduction of the growth rate there which de- 
lays the onset of turbulence. Keplerian flows are therefore 
expected to become turbulent for Reynolds numbers a few 
order of magnitudes higher than for plane subcritical flows 
(see: Longaretti 2002; Tevzadze et al. 2003). The possibility 
of an alternate route to turbulence gave new impetus to the 
research on the dynamics of astrophysical disks (Lominadze 
et al. 1988; Richard & Zahn 1999; Richard 2001; loannou & 
Kakouris 2001; Tagger 2001; Longaretti 2002; ChagclishviH 
et al. 2003; Tevzadze et al. 2003; Klahr & Bodenheimer 
2003; Yecko 2004; Afshordi et al. 2004 Umurhan & Regev 
2004; Umurhan & Shaviv 2005; Klahr 2004; Bodo et al. 
2005; Mukhopadhyay et al. 2005; Barraco & Marcus 2005; 
Johnson & Gammie 2005a, 2005b; Umurhan 2006). By 
adapting the progress of the hydrodynamic community to 
the disks flow, this research is promising for solving the 
disks' hydrodynamic turbulence problem. 

But it remains to be seen whether this route to turbu- 
lence actually applies to astrophysical disks. Compared to 
plane shear flows, these possess two additional properties: 
differential rotation and vertical stratification. Separate 
studies of these factors show that each exerts a stabiliz- 
ing effect on the flow: these include numerical calculation 
of the stability of unstratified flows by Shen et al. (2006), 
experiments on Keplerian rotation without stratification by 
Ji et al. (2006), estimates of the growth rates with strat- 
ification by Brandenburg & Dintrans (2006). However, it 
appears that the combined action of differential rotation 
and stratification introduces a new degree of freedom that 
may influence the flow stability and lead to turbulence at a 
high enough Reynolds number. Indeed, it has been shown 
that strato-rotational flows may exhibit global instability 
in bounded domains (Dubrulle et al. 2005). However, it is 
probable that local disk dynamics will also lead to hydro- 
dynamic turbulence. The study of the linear perturbations 
in strato-rotational flow in the local limit can be found in 
Tevzadze et al. (2003; hereafter T03); it is shown there that 
the combined action of rotation and stratiflcation generates 
an aperiodic vortex mode, which undergoes nonmodal tran- 
sient growth. We conjecture that this transient growth may 
be the main energy source for the turbulence in the bypass 
scenario. 

Although the importance of the transient exchange of 
energy between perturbations and mean flow has now been 
realized by most working in the field, only a few seem 
aware that another linear process may play also a crucial 
role, namely the linear coupling of modes, which allows for 
transient exchange of energy between them. As shown by 
Chagelishvih et al. (1997a,bj and Gogoberidze et al. (2004), 



the energy exchange between modes is inherent to shear 
fiows (as the transient exchange of energy between the mean 
flow and perturbations), and it determines in many respects 
the diversity of perturbation modes and, therefore, of non- 
linear processes. Once one fully realizes the role of nonlinear 
processes in the bypass concept discussed above, it becomes 
evident that the neglect of linear mode coupling may lead 
to an incorrect picture of nonlinear (and, consequently, tur- 
bulent) phenomena. 

There are signs of such linear mode coupling in the sim- 
ulations performed by Klahr (2004), Barraco & Marcus 
(2005), and Brandenburg & Dintrans (2006), but appar- 
ently they were not identified as such. Compressive waves 
arc present in the simulation by Johnson & Gammie 
(2005b), along with vortical perturbations, but their ori- 
gin is not recognized, namely linear mode coupling. On 
this mode coupling, attention is focused in T03 and Bodo 
et al. (2005). The latter paper studies the linear dynam- 
ics of an imposed two-dimensional pure vortex mode per- 
turbation, in a compressible Keplerian disk with constant 
mean pressure and density. (Two-dimensionality, i.e. the 
neglect of cross-disk variation, is only correct for pertur- 
bations with characteristic scales close to or larger than 
the vertical stratification scale of the disk.) Two modes - 
a vortex mode and a spiral-density wave mode - exist in 
the system, and they are strongly coupled. This investiga- 
tion points out the importance of mode coupling and the 
necessity of considering compressibility for dynamic pro- 
cesses with characteristic scales close to or larger than the 
disk thickness. 

In T03 we studied the linear dynamics of three- 
dimensional small-scale perturbations (with characteristic 
length scales much shorter than the disk thickness) in com- 
pressible, vertically (stably) stratified Keplerian disks. The 
first novelty presented in that paper is the occurrence of 
an interplay between the disk rotation and stratification. 
Separately, each of these factors is stabilizing. However, 
their interplay gives rise to a new vortex/aperiodic mode 
that is able to extract the basic flow energy transiently. The 
second novelty is the existence of a linear coupling of that 
vortex mode with spiral-density wave (SDW) modes, that 
makes SDWs valuable participants of the dynamical pro- 
cesses. Furthermore, we compared the linear dynamics of 
the small-scale perturbations in steady stratified disks with 
that of the perturbations in imboimdcd plane Couctte flow. 
We showed that just the linear coupling of the disk modes 
makes the dynamics in the disk and plane flows similar to 
each other (provided the Reynolds number of the disk flow 
is chosen about three orders of magnitudes higher than in 
the plane flow). This similarity suggests that the bypass 
concept can also be applied to disk flow and it further mo- 
tivates investigation in this direction. 

In this respect, the present paper is a sequel of T03: 
we focus again on the mathematical and physical aspects 
of mode coupling, while introducing a significant simpli- 
fication by neglecting the rotational-acoustic waves. This 
is justified by the fact that the characteristic timescale of 
these modes is much shorter than for the two other modes, 
when the characteristic lengthscale of the perturbations is 
much shorter than the disk thickness. Then the rotational- 
acoustic waves are not coupled to the vortex and SDW 
modes, and they play a negligible role in the slow, small- 
scale dynamics. That is why we can cut out the rotational- 
acoustic waves (i.e., the flow compressibility) without detri- 
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ment to the dynamic picture and confine ourselves to the 
Boussinesq approximation. Keeping just vortex and SDW 
modes, this approximation simplifies mathematical aspects 
of the problem and allows an advance in the analytical de- 
scription and comprehension of mode coupling. 

The paper is organized as follows. In Sect. 2 we in- 
troduce the physical approximations and the mathemati- 
cal formalism, and describe the linear strato-rotational bal- 
ance and the perturbation modes. In Sect. 3 we present the 
qualitative and quantitative analysis of the linear dynamics 
of perturbations. We summarize and discuss the results in 
Sect. 4. 



2. Disk model 

The dynamics of a rotating flow in a central gravity field is 
governed by the Navier Stokes equations, which are written 
here in the cylindrical coordinates: 



dVr 



IdP 



+ (V . V)K - ^ = -- „ „ . 
ot r p or or 

P = P{P,S), 
where 



(1) 

(2) 

(3) 

(4) 
(5) 

(6) 



is the component in convective derivative in cylindrical 
coordinate system. For the equilibrium state we choose 
the thin disk approximation with Vq = (0, Vo0,O), with 
= (0,0,f2(r)) and Vq^ = rn(r). Hence, we neglect self 
gravity and assume that the disk is rotationally supported. 
Vertical gravity acceleration and stratification scaleheight 
can be defined as follows: 



a$(r) 
dz 



= n'^{r)z « sgn{z)g 



1 dPoiz) 1 dpoiz) 



Po{z) dz po{z) dz 



-kn- 



(7) 
(8) 



We consider isothermal equilibrium and study small-scale 
perturbations employing local approximation for further 
analysis. 

2.1. The shearing sheet formalism in the Boussinesq 
approximation 

For the purpose of the linear analysis, we employ local co- 
rotating shearing sheet: 



x = r-rQ, y = ro{(f)-Q.ot), 

X y z 

— — — «r 1 

ro ro ro 



(9) 
(10) 



where we neglect the effect of global flow curvature and 
study the local influence of the differential rotation on the 
perturbations (see Goldreich & Linden-Bell 1965). We in- 
troduce the linear perturbations as follows: 



V = Vo + V, P = Po + P', p = po+p' . 



(11) 



For further simplification we employ the Boussinesq ap- 
proximation to neglect compressibility effects. We assume 
that kH = const and restrict analysis to the perturbations 
with shorter lengthscales than the vertical stratification 
scale of the disk fiow. Hence, the system of equations that 
governs the perturbation dynamics in the considered limit 
reduces to 



1 dP' 
" po ox 



+ 2Ax^^V^ + 2{no + A)V^ 



— 2Ax— 
dt dy 



T./ 1 dP' p' 
Po oz Po 



0, 

1 dP' 
Po dy 

= 0, 



= 0, 



d 



d 



- + 2Ax— ^ - knVi 
dt dyj Po 



dK^dv^^dV[ 

dx dy dz 



0. 



(12) 
(13) 
(14) 
(15) 
(16) 



Here for brevity po = Po(0), and only the z > case is 
considered. We use standard Oort's constants to quantify 
the radial velocity shear: 



A- 



dn 

dr 



-no, and A + B = -Q.o. (17) 



Following the standard procedure of nonmodal analysis, we 
employ spatial Fourier expansion of perturbations with spe- 
cific phase variation: 



P'(r,t) 

P'(r,i) 
V'(r,i) 



p(k,i) 
e(k,t) 
v(k,t) 



where 

kx{t) = fcx(O) - 2Akyt . 



exp{ika;{t)x + ikyy + ikzz), (18) 



(19) 



This leads to the linear system of ODE that governs the 
dynamics of SFH in the considered flow: 



dt 



Vx(t) = 2VloVy{t) - \kx{t) 



Pit) 



po 



2BVx{t)-iky^, 

po 



Vzit) = -ik, 



p{t) git) 



dt 

d g{t) 



Po 



9- 



Po 



kHVz{t), 



dt Po 

{t)Vx{t)+k 

y^y ~^ ^z'^z 



(20) 
(21) 
(22) 

(23) 

(24) 
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Introducing new notations for pressure and density 



Po 



Po 



(25) 



we obtain an ODE system that does not explicitly contain 
complex coefficients. The system is characterized by a tem- 
poral invariant that corresponds to the linear perturbation 
of the potential vorticity: 

k 

W = kr,{t)vy ~ kyV^ + 2B--^D = const. (26) 

kn 

With straightforward simplifications 

k 2Qn k 
P = -^-j^^ + "fc^ ik4t)vy - kyv^) - 4A^w, , (27) 

we reduce the governing equations to the following ODE 
system that is third order in time: 



-ri^xW 



(28) 



2(.o + 2A)M^.. + 2.„^.,+,M|.^ 



d_ 

At 



17 "a W 



(29) 



op fcx(0 + I f) ;| \ „, 90 kx{t)ky kykz 

k^{t) + ^^W) ] ""^-^^'^W '^W) 



with 

Vz{t) = {kx{t)Vx{t) + kyVyit)) 



(30) 



(31) 



These equations will be employed in the numerical calcu- 
lations. The spectral energy of the perturbation SFH can 
be derived from the sum of thermal and thermobaric ener- 
gies (see Eckart 1960), which in the considered Boussinesq 
approximation will read as 



Ekit) = f (ulit) + ulit) + ul{t) + ^D\t) 



(32) 



is the Brunt- Vaisala frequency that is constant in our for- 
malism. The solutions to this dispersion equation consist of 
two modes that are generally decoupled in the zero shear 
limit: 

i) A spiral-density wave (SDW) - or inertia-gravity 
wave - that describes the oscillations due to the Coriolis 
force and to the vertical buoyancy; 

ii) A stationary vortex mode with zero frequency. This 
stationary solution describes aperiodic perturbations to the 
vertical strato-rotational balance that can be directly de- 
rived from Eqs. (|12m6p : 



2r2o— curiy 
oz 



92 



(35) 



where curly' is the vertical component of the vorticity: 



curiy = — ^ 
ox 



dy 



(36) 



We see that stationary vorticity perturbations in the disk 
plane (lefthand side of the Eq. [55]) can be sustained by den- 
sity perturbations (righthand side of the Eq. [35]), provided 
that both rotation (fio) Etnd vertical stratification (g) are 
present. It is the stationary vortex mode, which appears 
due to the strato-rotational balance, that becomes dynam- 
ically active in the shear case (for A ^ 0) and is then able 
to extract the mean flow energy. 

3. Qualitative analysis of the linear dynamics: 
transient growth and mode coupling 

The numerical study of SFHs dynamics is governed by 
Eqs. (|28ll3ip . However, to clarify the physical nature of the 
perturbations and their linear dynamics, it is advisable to 
rewrite them in another form. Introducing a new variable 



(37) 



from Eqs. (27),(29)-(31) one obtains the second-order in- 
homogeneous differential equation 



di2 



(38) 



with the frequency 



2.2. Rigid rotation limit: strato-rotational balance and 
perturbation modes 

In the shearless limit (A = 0, B = — fi), the operators 
that intervene in the system above become normal, and it 
is easy to derive its dispersion relation. Its roots will help 
us to understand the mode behavior in the Keplerian flow, 
although the results are not applicable as such. Taking the 
Fourier expansion of the perturbations in time oc exp{iLut), 
we obtain 



OJ U) 



where 



/r2 







N' = gkH 



9 dpo{z) 
Pq{z) dz 



(33) 



(34) 



(t) = 7V2^-K4f^o(f^o + A)-f - 



and the coupling function 



(39) 



8^f7o-4A2^ + m2^^, 



(40) 



Equation (|38p is essentially the same as the main dy- 
namical equation of Chagelishvili et al. (1997), in which 
one form of mode linear coupling in shear flows has been 
described, namely the generation of wave SFHs by related 
SFHs of vortex mode perturbations. Thus one can borrow 
the qualitative analysis from that paper. As a result, Eq. 
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Fig. 1. The evolution of a single SFH when initial 
conditions correspond to the tightly leading pure vor- 
tex mode. Velocity components (vx{t)kH I^q, Vy{t)kH /^o, 
Vzit)kH /^o,) and normalized energy of perturbation SFH 
{Ek{t)/E{0)) are shown. Here kx{0)/ky = -30, ky = k,_ = 
lOfcff. 
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Fig. 2. Same as Fig. 1, but when the initial conditions for 
the single SFH correspond to a mix of vortex and SDW 
modes. Note that the differences are rather minor, since 
the vortex mode dominates. 



([55]) describes the dynamics of SFHs of two different modes 
of perturbations: 

(a) the spiral-density wave mode (SDW), which is de- 
scribed by the general solution of the corresponding homo- 
geneous equation. Consequently, this wave mode has zero 
potential vorticity (i.e. W — 0, cf. Eq. [^S]). Note that the 
frequency uj^{t) of the wave SFHs is time-dependent. 

(b) the aperiodic vortex mode, that originates from the 
equation inhomogeneity {f^{t)W) and is associated with 
the particular solution of the inhomogeneous equation. The 
amplitude of the vortex mode is proportional to fcf,{t)W 
and tends to zero with 0. In other words, the vortex 
mode has nonzero potential vorticity and at the same time 
acquires nonzero density perturbation. 

Any perturbation at any moment can be decomposed in the 
sum of a wave and a vortex mode. 

The character of the dynamics depends on which mode 
SFH is initially imposed in Eq. ([55]) : pure wave, pure vor- 
tex, or a mixture of the two. Note that the dynamics of 
tight SFHs {\kx{0)/ky\ » 1) is adiabatic, whereas that of 
open SFHs {\kx{0)/ky\ ^ 1) is non- adiabatic. Thus a SFH, 
which initially is tightly leading {kx{t)/ky -C —1), becomes 
open in due course (with non-adiabatic dynamics) , and will 
finally be tightly trailing {kx{t)/ky ^ 1). 



To impose initially pure vortex mode perturbations, we 
use approximate solutions that can be derived from Eq. 

*o(i) = . (41) 

This solution describes a vortex mode when the temporal 
variation of the vortical perturbations is much slower than 
the SDW oscillations. Hence, we start our numerical inte- 
gration of Eqs. (|28ll3ip at wavenumbers where the above 
condition is amply satisfied, i.e., at \kx{Q)/ky \ 3> 1. 

Figure [T] displays the evolution of a dynamically active 
{ky — kz ) single SFH when the initial conditions correspond 
to a tightly leading pure vortex mode {kx{0)/ky = —30). 
We see that in the leading phase only that vortex mode 
is present. Oscillations (hence wave SFH) appear in the 
trailing phase after time t — t* — 20/^0 when kx{t) = 0. 
The SDW oscillates mainly in the vertical direction; thus 
oscillations are the most pronounced in Vz- The graph of 
the density D is similar to that of Vz and so is not pre- 
sented here. The energy graph shows the transient growth 
in the leading phase. In the trailing phase, the generated 
wave SFH keeps the energy (since SDW do not exchange 
the energy with the mean flow). This course of events - the 
transient growth and subsequent coupling - makes the en- 
ergy dynamics similar to the 3D unbounded plane Couette 
flow case. Note that the tighter the initial leading SFH is, 
i.e. the higher the ratio \kx{0)/ky\, the stronger the tran- 
sient growth. The possible maximum value of \kx{0)/ky \ is 
determined by the Reynolds number (Re). Consequently, 
the amplification factor is determined by Re and it can 
become huge since in Keplerian disks Re reaches literally 
astronomical values {Re > 10^°) (see Sect. 3 of T03 for 
details) . 

Imposing a leading SFH {kx{0)/ky < 0) of pure wave 
mode, we do not detect any notable transient growth in the 
wave SFH, or any generation of the vortex mode: since the 
potential vorticity is time invariant (see Eq. I26p , the wave 
SFH, having zero potential vorticity, is not able to gener- 
ate the vortex SFH, which has non-zero potential vorticity. 
Hence, waves with zero potential vorticity alone, which were 
also recently described by Goodman and Balbus (2001), are 
not expected to support energetically hydrodynamic turbu- 
lence in Keplerian flows. 

Figure [2] shows the evolution of another dynamically 
active single SFH with ky — kz, again with kx{0)/ky = 
—30, but where the initial conditions correspond to a mix of 
vortex and SDW modes. The initial admixture of the wave 
SFH does not change the dynamic picture qualitatively, 
due to the already mentioned reason, that SDW does not 
exchange energy with the mean flow. The energy graph is 
similar to Fig. 7 of Brandenburg & Dintrans (2006), though 
the linear mode coupling has not been identified there. 

Let us concentrate on the linear dynamics when we ini- 
tially insert a tightly leading SFH of vortex mode pertur- 
bation. The dynamics is then the most complex and the 
energy exchange the most efficient, due both to the tran- 
sient growth of the SFH and to mode coupling with the 
SDW. The intensity of these processes depends on values 
of kx{0)/ky and kz/ky. 

The ratio Ek{t*)/E{0) gives a good estimate of the non- 
modal growth of vortex SFHs. Here, t* is the time when 
kx{t*) = 0. The spectral behavior (specifically, the depen- 
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Kz/Ky Kz/Ky 



Fig. 3. Left: energy amplification through transient growth 
of vortex mode perturbations Max{E {t) / E (0)) vs. k^/ky. 
Right: mode coupling efficiency measured by the generated 
wave energy E{wave)/ E{0) vs. kz/ky. Here kj;{0)/ky — 
—30, ky = 10 kn and the maximum efficiency of mode cou- 
pling is achieved at kz/ky = 1.25. 



dence on kz/ky) of the transient amplification rate calcu- 
lated by Eqs. (|28ll3ip is shown in the left panel of Fig. El A 
constant value of kx{0)/ky = —30 is chosen to simulate the 
energy growth of vortex perturbations at constant Reynolds 
number. One can see that the transient amplification is ef- 
ficient at low values of kz, while it is strongly reduced for 
SFHs with kz :^ ky. 

The right panel of Fig. [3] shows the mode coupling effi- 
ciency, measured in units Ek{wave)/ E(0), vs kz/ky, where 
Ek (wave) is the final energy of the wave SFH (Ek (wave) — 
E{tf^n) - E{0), tfin = 2t* = 4Gf7o). That efficiency has a 
pronounced maximum at kz — ky. The energy of the am- 
plified vortex SFH is transformed into the wave SHF due 
to mode coupling, so at later times, the main carrier of the 
perturbation energy is the SDW (since vortex decay in the 
trailing phase, while waves evolve keeping nearly constant 
energy). Therefore, SFHs with kz ^ ky are dynamically 
active and should play an important role in the turbulent 
process. These wavenumbers were chosen in Figs. [l]l2] to 
show the dynamical picture of these SFHs and illustrate the 
significance of the transient growth and subsequent mode 
coupling. 

In Fig. |4]we have chosen the phase space to display the 
evolution of the vertical and total velocities. The initial con- 
ditions correspond to the pure vortex mode (top row, A.l 
and A.2) and a mix of vortex and wave modes (bottom row, 
B.l and B.2) presented in Figs. [3]and|31 The initial point 
is marked by a triangle, and the point at which k^it) = 
is marked by a circle. The segments from triangle to cir- 
cle in the top row graphs describe the transient growth of 
the leading phase of the vortex mode. In the trailing phase, 
after the circle marker, one can see some subsiding of the 
vortex harmonic and generation of the wave (due to mode 
coupling). The latter is indicated by the cyclic phase trajec- 
tories. The persistence of cyclic trajectories illustrates the 
constancy of the wave SFH energy: the wave is the main 
carrier of the perturbation energy at later times. The bot- 
tom graphs show that the mix of initial vortex and wave 
harmonics evolves much in the same manner. Differences 
are seen initially, before the vortex has undergone tran- 
sient growth. The graph B.l shows the transient growth 
and mode coupling in phase space in terms of the vertical 
velocity. Note how the initial "low level" cyclic phase trajec- 
tories are transformed into "high level" cyclic trajectories 
by transient growth of the vortex admixture and the subse- 
quent generation of the wave harmonic by the vortex mode. 




Fig. 4. The vertical velocity (panels A.l and B.l) and total 
velocity amplitude (panels A.2 and B.2) of perturbations 
in phase space (v, dv/dt). Top row shows the evolution 
of perturbations when initial conditions correspond to the 
vortex mode (panels A.l and A.2), while a mix of vortex 
and wave modes is shown in the bottom row (panels B.l 
and B.2). A triangle marks the initial state, while the circle 
corresponds to the time when k,j;(t) = 0. 

It illustrates how transient growth raises a perturbation to 
a higher energy state. 

The vortex mode undergoes the amplification or damp- 
ing depending on its phase (leading or trailing) in differ- 
entially rotating fiows. From this, it has been speculated 
(see Johnson & Gammie 2005) that any realistic ensem- 
ble of these modes covering a wide area in wave-number 
space will not exhibit cumulative energy growth. However, 
as we have seen in Figs.[T][2]total energy of the SFH exhibits 
asymmetry in growth and decay, due to the mode coupling. 
Hence, any ensemble of initial perturbations that include 
vortex mode perturbations and cover the region kz ^ ky 
will undergo overall energy growth in time. 

4. Discussion 

Our motivation in studying mode coupling in the presence 
of vertical stratification is quite clear: in unstratified rotat- 
ing flows the strato-rotational balance is absent (see Eq. 
I35p . hence there will be no vortex mode, whose role is so 
powerful in extracting the shear flow energy. Moreover, in 
stratified rotating fiows mode coupling causes the gener- 
ation of SDW, which are able to conserve the extracted 
energy. Thus one expects astrophysical disks, which are ver- 
tically stratified, to demonstrate intrinsically different be- 
havior compared to unstratified rotating flows. Therefore it 
is not possible to draw conclusions about the stability of as- 
trophysical disks by considering unstratified rotating flows, 
as has been done often in the past and again in recent in- 
vestigations: analytical (Balbus 2006), numerical (Shen et 
al. 2006), and experimental (Ji et al. 2006). 

Our analysis shows that, in the local limit, spectrally 
stable stratified Keplerian disks allow for two modes of 
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perturbations - vortex and SDW - that are jointly able 
to extract the background flow energy and determine the 
disk dynamical activity in the small-scale range. These 
modes are linearly coupled due to the non-normality of 
Keplerian/shear flow. The coupling is asymmetric: the vor- 
tex mode is able to generate the related SDW, but the 
inverse is not true. This mode coupling is transient (like 
the energy exchange between perturbations and the basic 
flow): the SFH of the vortex mode generates the wave SFHs 
during the brief time interval where it switches from leading 
to trailing, thus rendering the dynamics non-adiabatic. 

At first sight, the mode coupling described here seems 
similar to the phenomenon of geostrophic adjustment, 
since both lead to the generation of the spiral density 
waves. However, these processes are intrinsically different. 
Geostrophic adjustment is an initial value problem that de- 
scribes the transition from an initially unbalanced state 
to that of geostrophic balance. In the general case, part 
of the initial conditions containing no potential vorticity 
(the wave component) will radiate away as inertial-gravity 
waves (see, e.g., Pedlosky 2003) during the geostrophic ad- 
justment, whereas our initial conditions do not include zero 
potential vorticity corresponding to the wave component. 
Hence, we have eliminated the wave generation due to the 
process of the initial geostropic adjustment. The waves gen- 
erated in our case stem from linear mode coupling induced 
by the velocity shear. Moreover, geostrophic adjustment is 
mainly a nonlinear process that leads to equilibrium. In 
contrast, wave generation due to mode coupling is a linear 
process and does not describe the relaxation of the system, 
but hopefully the opposite: it promotes the transition to 
turbulence. 

The linear dynamics of each leading SFH of vortex mode 
proceeds in the following way: initially, the SFH extracts 
energy from the basic flow and it grows. At the same time 
kx{t)/ky — > —0. Becoming trailing, the vortex SFH gener- 
ates the related SFH of SDW. In what follows, while tilting 
(i.e., increasing kx{t)/ky), the wave SFH keeps the energy, 
whereas the vortex SFH gives its energy back to the basic 
flow, so the energy gained by the leading vortex SFH is con- 
served by the SDW. This course of events - the transient 
growth plus coupling - is strongly pronounced for SFHs 
with kz/ky ~ 1 and makes their energy dynamics similar 
to that of a 3D unbounded plane Couette flow. 

We are aware that there are differences between them: 
in the Couette case, only the vortex mode participates in 
the dynamical process, whereas in the disk case, this role 
is played (in the local limit) by the symbiosis between vor- 
tex and SDW perturbations. However, given the similari- 
ties that we have discussed, we conjecture that the bypass 
concept, which has been developed for the transition to tur- 
bulence in the Couette flow, is also applicable to rotating 
stratified disks. One question remains, of course: will the 
nonlinear interactions provide positive feedback that is effi- 
cient enough for that bypass transition? To answer it, more 
numerical simulations are needed, which must include the 
vertical stratification and grasp the mode coupling. 
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